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Abstract 

We study effects of the medium polarization on superfluidity in symmetric nuclear matter in a 
relativistic formalism. An effect of the liquid-gas instability is emphasized. We examine two types 
of decomposition of the nucleon propagator; the standard Feynman-density and the particle-hole- 
antiparticle ones. In both cases, the medium polarization effect is determined by a characteristic 
cancellation among the a, the longitudinal uj, and the a-uj mixed polarizations. The instability 
leads to increase of pairing gap. Around the saturation density that is free from the instability the 
medium polarization enhances pairing gap in the former case whereas reduces in the latter. At the 
lowest density that is also free from the instability the gap increases in both cases. 
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I. INTRODUCTION 



Superfluidity in nuclear matter has long been studied mainly in pure neutron matter 
from a view point of neutron-star physics such as cooling rates and glitch phenomena Q]. In 
addition, that in nuclear matter with finite ZjN ratio is also becoming of interest as basic 
information for the structure theory of finite nuclei, since recent development of Rl-beam 
experiments makes it possible to study ~ Z medium-heavy nuclei and neutron-rich light 
nuclei. 

At present, there are two ways to describe the fundamental properties such as the satu- 
ration property of the finite-density nuclear many-body system; the non-relativistic and the 
relativistic models. They are understood as describing observed properties almost equally. 
Among them, we here adopt the latter. The origin of quantum hadrodynamics (QHD) can 
be traced back to Duerr's relativistic nuclear model ^ that reformulated a non-relativistic 
field theoretical model of Johnson and Teller Q|. Since Chin and Walecka succeeded in 
reproducing the saturation property of symmetric nuclear matter within the mean-field ap- 
proximation , QHD has not only been evolving beyond the mean-field approximation 
as a many-body theory but also been enlarging its objects as a nuclear structure model such 
as infinite matter — spherical deformed — > rotating nuclei These successes indicate 
that the particle-hole channel interaction in QHD is realistic. In contrast, relativistic nuclear 
structure calculations with pairing done so far have been using particle-particle channel in- 
teractions borrowed from non-relativistic models and therefore the particle-particle channel 
in QHD has not been studied fully even in infinite matter. Aside from practical successes, 
this situation is unsatisfactory theoretically. Therefore, in this paper, we present an effort 
to derive an in-medium particle-particle interaction that is consistent with the relativistic 
mean field (RMF) although only infinite matter can be discussed at the present stage. 

Up to now, there have been a lot of non-relativistic studies of pairing in nuclear matter. 
As for the particle-particle interaction entering into the gap equation, many authors adopted 
bare interactions whereas others adopted renormalized ones such as G-matrices. Although, 
in the medium, renormalized interactions would be used intuitively, following reasons support 
the use of bare interactions: (1) The Green's function formalism leads to the sum of the 
irreducible diagrams 0, 1^, and its lowest order is the bare interaction. (2) The gap equation 
itself implies the short-range correlation U, In general, medium renormalizations 
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are expected to enhance the gap by reducing the short-range repulsion. An interesting 
exception is the Gogny force. This is known to reproduce the pairing properties given by 
bare interactions at least at low densities Anyway, as a next step, polarization diagrams 
should be considered. In the non-relativisitc framework, a lot of works have been done to 



study medium polarization effects on superfluidity in neutron matter 



IS 
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Q, 15, 16, 0, Q, 



. All of them concluded that the medium polarization reduces pairing gap 
significantly. A recent ab initio calculation [23^, however, claimed that the polarization 
effect is weak. Studies of polarization effects on symmetric nuclear matter have just begun 
recently. Reference j^] reported that the gap increases substantially. At the low density 
limit, Ref. j3] discussed that in the Fermi systems with 4 species the medium polarization 
enhances the gap. 

Symmetric matter was also studied in the relativistic framework 0, 26 1. Both of these 



references reported that the vacuum polarization reduces the gap while the latter reported 



that the medium polarization enhances the gap. References [2f 
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discussed that the 



increase of the gap is related to the existence of an instability. This is known as the liquid- 



gas instability [2? 



Thus, in this paper, we investigate the liquid-gas instability and its effects on superfluidity 
in symmetric nuclear matter. Results on pure neutron matter are also mentioned briefly. 
Preliminary results were reported in Ref. j3l|. 



II. LIQUID-GAS INSTABILITY IN RELATIVISTIC RANDOM PHASE AP- 
PROXIMATION 

We begin with the ordinary a-uj model Lagrangian density, 

C = ^{i-f^d" - M)^ 

^^lu = df,uJu - dyujf, . (1) 

Here ■?/', a, and uo are the nucleon, the a meson, and the uj meson fields, respectively, M, 
mo-, and are their masses, and and (7^ are the nucleon- meson coupling constants. 
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The relativistic mean field approximation is carried out by replacing the meson fields in the 
coupled equations of motion by their expectation values as 



a ^ {a) = ao, 

i {<^ti) = Sfio^^o ■ (2) 



Then the nucleon effective mass (Dirac mass) equation is given by 

M* = M- g^ao 

where the isospin factor A = 2 and 1 indicate symmetric nuclear matter and pure neutron 
matter, respectively. The Fermi momentum is related to the baryon density as 



PB = ■ (4) 



Properties of normal fiuid, zero temperature matter of a given density is completely de- 
termined by the above effective mass equation Q. The so-called saturation curve or the 
equation of state is given by the binding energy per nucleon, E/A — M = S/pb — M, as a 
function of pe or kp. This immediately gives pressure. 



2 d f 8\ 



dpB \pbJ 

Thermodynamic stability of the matter in the liquid phase is stated as dP/dp-Q > 0. Since 
the ratio of the sound velocity Cg to the light velocity c is given by 



(cg/c)^ < in FigH indicates the existence of a mechanical instability to the gas phase. In 
this pci^er we adopt M = 939 MeV, = 550 MeV, = 783 MeV, gl = 91.64, and g^ = 

136.2 y. 

Quantum mechanically, stability of a state is determined by the second variation of the 
energy with respect to the fields 28|. This is equivalent to the random phase approximation 
(RPA). The RPA in the present model in which the nucleon- nucleon interaction is mediated 
by mesons is formulated by calculating the meson propagators that couple to the particle- 



hole and particle-antiparticle polarizations 



30|. The Dyson equation that determines 




FIG. 1: Squared ratios of the sound velocity to the light velocity as a function of the Fermi 
momentum. 

the RPA propagator D is given by 

D = D^ + D^YiD (7) 
(pictorially given in Fig|2)), where the lowest order propagator, 

= I 1 I . (8) 



and the polarization insertion, 



are given by 5x5 matrices. Their components are specified as 







n= : , (9) 



q2 _ _|_ 



and 



KM = (g,. - ^) D^iQ) , D^iq) = ^ , (10) 

KM = -^al J -0j-Mi,G{khMk + q)] , 

K^q)=W.9.j (0iTr[7.G(fc)G(fc + g)] . (11) 

Here G{k) stands for the nucleon propagator and the Tr symbol includes isospin. IlJ^ stands 
for the matter induced a-uj mixed polarization; a a excites a particle-hole pair and then it 
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decays into an uj and vice versa. Since Eq.(|7j) is formally solved as 

1 



D= ^ \^ Do, (12) 



zeros of the dielectric function, 



determines collective excitations. 



e = det (1 - DoU) (13) 



FIG. 2: Feynman diagram representing the RPA meson propagator. 

For the present purpose, investigating instability, only the real parts at zero energy trans- 
fer are necessary. Therefore we set 

-1 



|q|z _|_ 

^o'' = I |,i , , (14) 
|q|^ + 

and the second term in -D^^ drops because of the baryon number conservation. This con- 
servation law also restricts non- vanishing components among 11; if we choose the coordinate 
system as g = (g°,0,0, |q|), only = U^q - and = = among U^^ and 
n" = among 11^ survive. Note that energy transfer is set to zero, that is the instan- 
taneous approximation. After some permutations, 1 — DqU becomes block diagonal, and 
consequently the dielectric function reduces to 



eL = (1 - D^,U^){1 - Z^o'n^) - D',D',{Uy , 

CT = 1 + D^U^ . (15) 

The transverse dielectric function is always positive in the density region in which su- 
perfluidity is realized, whereas the longitudinal one eL becomes negative at intermediate 
densities; this represents the liquid-gas instability. 

In this work, we concentrate on the particle-hole polarization. We examine two ways to 
extract the particle-hole polarization; (1) the standard Feynman- density (FD) decomposition 
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of G{k) and (2) the particle-hole-antiparticle (pha) decomposition. From its definition, tlie 
nucleon propagator in the medium is given by 

i-e{kF-\k\) e{kF-\k\) 



- i^^K^' + M* 



2E*{k)v'^ '\k^ - E*{k) + it k^ - E*{k) - ie 

1 



where 



k^ + E*{k) -ie\ 
Gp{k) + G^{k) + G,{k) , 



K={E*ik),k), k={-E*ik),k), 



(16) 



E*{k) = Vk2 + M*2. (17) 

The first, the second, and the third terms represent the propagator of particle, hole, and 
antiparticle, respectively. By sorting them with respect to the Heaviside function, another 
form, 



Gik) ^ i,,k^ + M-n (^,3^^ + ^^(^° - E^kmk. - iki) 

= GM + Gu{k) , 



(18) 



is obtained. The first and the second terms are called the Feynman and the density parts, 
respectively. The former consists of the antiparticle propagation and a part of the particle 
propagation, while the latter consists of the hole propagation and the other part of the 
particle propagation. In the standard FD decomposition, the density dependent GpGo + 
GbGf part in Eq. (|TT| is regarded as the particle-hole polarization. Note that the GbGb 
part is pure imaginary and therefore is not necessary for the present purpose. In the pha 



e-hole polarization 
The concrete 



decomposition, the GpG^ + Gi^Gp part in Eq. (fTT| represents the partic. 
and this directly corresponds to that in non-relativistic calculations 
expressions of 11 are given in Appendix A. 

Figure El shows the longitudinal dielectric function eL of the FD and the pha cases for 
symmetric nuclear matter. At low momentum transfers, eL becomes negative in both cases. 
The density range of the instability coincides very well with Fig. ^ in the FD case. This is 
consistent with Ref. j3] in which the density range of the instability hardly changes even 
after inclusion of the vacuum polarization. In the pha case the instability region shrinks. 
These results indicate that the liquid-gas instabihty occurred in a wide density range affects 





FIG. 3: (Color online) Longitudinal dielectric functions for symmetric nuclear matter as a function 
of the Fermi momentum and the momentum transfer, (a) the FD case, (b) the pha case. 

the pairing properties that are calculated by using the RPA meson propagator that exhibits 
instability. 

III. MEDIUM POLARIZATION EFFECTS ON SUPERFLUIDITY 



The concrete form of the RPA meson propagator is given by 



D 



1 



1-Don 



Do 



V 









EL 







ex 









^0 





\ 







This gives the particle-particle channel interaction that determines superfluidity as 



ga 



a,b 



1 : a = -1 
7m ■ a = 



(19) 



(20) 



The antisymmetrized matrix element of this interaction for the ^5*0 pairing channel is given 
by 

v{p, k) = (ps', ps'|VRPA|ks, ks) - (ps', ps'| VRPA|ks, ks), (21) 
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here the argument |q| of D in Vrpa is specified as q = p— k, and tildes represent time reversal. 
Since, by ignoring the coupling to the negative-en ergy states, the relativistic Numbu-Gor'kov 



formalism reduces to the usual gap equation 



after an integration with respect to the angle between p and k. The effective mass equation 
(jni) is slightly modified when superfiuidity sets in as 



M* = M-^— vjk'dk . 

1 /. Ek — Ekp 



Here E^ = E*{k) + Quji^^o)- Using these equations we calculate superfiuidity in nuclear 
matter assuming that it is in the liquid phase even at the density range in which the liquid- 
gas instability would occur, as usual. We set the upper bound of the integrations as 20 
fm~^. 

Here, some discussion about our choice of the interaction is in order. In the lowest order 
(the tree level; D Dq), Vrpa reduces to the one boson exchange (OBE) interaction Vobe 
given by the RMF vertices (see Figip. It has been well known that this OBE interaction 
gives unphysically large pairing gap jSJ]. Its reason can be traced back to the fact that the 
RMF vertices were tuned only below the Fermi momentum. In order to discuss the higher 
order (polarization) effects this OBE result must be improved beforehand. We accomplish 
this by introducing a form factor that modulates the high momentum part of the interac- 
tion smoothly as follows. Note that the authors of Ref. [3| adopted sudden momentum 



cutoffs so as to reproduce the virtual state in the T matrix. We argued in Ref. |36| that the 

n 

sudden cutoff distorts the shape of the short range pair wave function. See also Ref. j37|] for 
the sudden momentum cutoff that reproduces the results of a bare interaction in the pair- 
ing calculation. From a general argument, the lowest order in the particle-particle channel 
interaction that gives pairing should be a bare interaction and the particle-hole channel in- 
teraction that determines the medium polarization should be an in-medium one in principle. 
In the present investigation, however, we calculate both the tree and the bubble contribu- 
tions on the same footing adopting an interaction of in-medium nature, which reproduces 
the results of a bare interaction in the tree level. That is, we regard, in a sense, the present 
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RMF-OBE interaction with a form factor resembles the Gogny force in the non-relativistic 
pairing calculations. 

In order to modulate the high momentum interaction that enters into the pairing calcu- 
lation, we introduce a form factor, 

at each vertex. This does not affect the mean field (Hartree) part with the momentum 
transfer q = 0. The parameter A is determined so as to minimize the difference in the pairing 
properties from the results of the RMF+Bonn calculation, that is, a hybrid calculation 
performed by adopting single particle states from the RMF model and the Bonn potential 
as the pairing interaction. Here we adopt the Bonn-B potential because this has a moderate 
property among the available (charge-independent) versions A, B, and C [38]. The pair wave 
function. 



2Eqp(fc) 



E^^{k) = ^{Ek-E,^y + A%k), (25) 

is related to the gap at the Fermi surface, 

A{kF) = -/ vikF,k)(l){k)k^dk, (26) 

4vr^ Jo 

and its derivative determines the coherence length, 

which measures the spatial size of the Cooper pairs. These expressions indicate that A(A;f) 
and ^ carry independent information, (f) and respectively, in strongly-coupled systems, 
whereas they are intimately related to each other in weakly-coupled ones. Therefore we 
search for A that minimizes 

2 ^ \ 2" 



2 ^ \ - I I A(^f)rmf - A(A;f) Bonn \ , / ^RMF ^Bonn \ i ,r,a\ 

^ "^irlV ^(^^ ) V ^Bonn ) /' ^'^^ 

with N = 11 {kp = 0.2, 0.3,. ..,1.2 fm"^). The obtained value is A = 7.26 fm'^ 0. For 
consistency we included the form factor also in the polarization diagrams but it almost would 
not affect them because the polarization affects only low momenta (see following figures). 
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FIG. 4: Cross sections of the dielectric functions for the FD case of symmetric nuclear matter as 
a function of the momentum transfer, (a) kp = 0.8 fm^-*^, (b) = 1.4 fm^-*^. 

Equation (|TTI|l indicates that Vrpa becomes ill-defined when the liquid-gas instability oc- 
curs. This means that superfluidity in the liquid-gas coexistent phase rather than that in the 
liquid phase should be considered. However, in the present calculation, we consider the case 
in which the system stays in the liquid phase also in the instability region (0.3 fm~^ < k-p < 
1.3 fm~^ in the FD case) as usual. Therefore in that region only qualitative discussion is 

n 

possible. Reference ^21 also mentioned the existence of the instability. Note that discussion 
is quantitative in high and low density regions that are free from the instability. Figure ID 
shows the cross sections of the longitudinal and the transverse dielectric functions. Fig- 
ureEfa) is for kp = 0.8 fm~^ and (b) is for kp = 1.4 fm~^. The former indicates that a low 
momentum cutoff is necessary to regularize the calculation for the instability region. Thus, 
we introduce ecut that operates to cut |q| for which eL/T(0, |q|) < ecut- Since other param- 
eters were determined at the saturation density, we chose ecut = 0.65 that maintains the 
full variation of eL and ex around this density (see FiglD^b)). This also serves to make the 
/cp-dependence of A{kp) smooth at the boundary of the instability region. Figure El reports 
the dependence of A{kp = 0.8 fm^^) on Ccut- Its dependence is moderate around the chosen 
value. Figure El compares A{kp). This figure shows that the medium polarization increases 
the gap at all densities in the FD case. This is common to the previous calculation 
in which the form factor to modulate the high momentum interaction was not introduced. 
Batata e. al. Q w.o weaWd effect ot .he bubble cba,.. b. tbe. . pa_ 
obtained a similar result. The contents of the polarization interaction, Vrpa — Vqbe; are 
decomposed in Fig.[7| (The /c-dependence in Fig. [Tfa) is a result of ecut for q.) These figures 
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represent a characteristic feature that the a polarization and the longitudinal uj polarization 
give strong attractions whereas the a-uj mixed polarization gives a strong repulsion and 
they strongly cancel each other. Remaining tiny attraction leads to the increase of A{k-p). 
From this viewpoint, the neglect of the a-u mixed polarization in Ref. would cause 
an imbalance. The transverse uo polarization that represents the spin density fluctuation 
is slightly repulsive, but the repulsion in the momentum region in which A(A;) < (that 
is predominantly determined by Vobe) also increases ^{k-p) because of the structure of the 
gap equation (j221)- 




I • ■ • • 1 

0.2 0.4 0.6 0.8 1 



FIG. 5: Cutoff parameter dependence of A(A;f) at kp = 0.8 fm 
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FIG. 6: (Color online) Pairing gaps at the Fermi surface for symmetric nuclear matter as a function 
of the Fermi momentum. 

Figure IHl also reports the pha case. The results of the FD and the pha cases are very 
close to each other at low densities. This is because the particle propagation contained Gd 
is small. However, their difference grows as density increases; at high k^ the polarization 
reduces A(/cf) in contrast to the FD case. This is brought about by the behavior that the 
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FIG. 7: (Color online) Decomposition of the polarization interaction for the FD case of symmetric 
nuclear matter as a function of the momentum, (a) kp = 0.8 fm~^, (b) kp = 1.4 fm~-^. 

total polarization interaction becomes repulsive (Fig. IH]). 




FIG. 8: (Color online) The same as Fig. d but for the pha case. 

Finally we briefly mention the pure neutron matter case. It is subtle whether liquid-gas 
instability occurs in pure neutron matter or not. In the present calculations it occurs in the 
FD case (Fig. IHl^a)). In the pha case it does not occur but eL strongly decreases at medium 
kp (Fig. in^b)). Consequently the behavior of A{kp) is similar to the symmetric matter case 
(Fig.Cni). 

IV. DISCUSSION AND CONCLUSION 



In symmetric nuclear matter, the RPA leads to liquid-gas instability at moderate densities 
even at zero temperature. In the present study, we concentrated on the medium polarization, 
but Ref. [3] showed that liquid-gas instability survives even after inclusion of the vacuum 
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FIG. 9: (Color online) The same as Fig. but for pure neutron matter. 
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FIG. 10: (Color online) The same as Fig. but for pure neutron matter. 



polarization. This indicates that the medium polarization effects on superfluidity should 
be considered in the liquid-gas coexistent phase rather than in the liquid phase. In the 
present investigation, however, we have studied superfluidity in the liquid phase, as usual. 
Therefore only qualitative discussion is allowed for the density region in which the instability 
occurs; the gap increases from the OBE result. In other regions free from the instability, 
a quantitative discussion is possible; the gap increases in the FD decomposition whereas 
decreases in the pha decomposition at high densities. At low densities the gap increases 
in both cases. These results are brought about by the characteristic cancellation between 
the attraction from the a polarization and the longitudinal lo polarization and the repulsion 
from the o-lo mixed polarization. The result for the FD case is consistent with that of 



another relativistic study in Ref. j3l and of a non-relativistic study in Ref. 27 1. Both cases 
are consistent with Ref. It is known that the coupling to surface vibrations enhances 
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pairing in finite nuclei [39|, |40|, |41[ . This tendency is in the same direction as the FD case 
in the present calculation. But, since the polarization in the matter case is density modes, 
their mutual correspondence should be clarified. 

It is a subtle problem whether liquid-gas instability occurs in pure neutron matter or 
not. In the present calculation, it occurs in the FD case. Although it does not occur 
in the pha case, the longitudinal dielectric function decreases strongly. Consequently the 
results for the gap are similar to the symmetric matter case. Note that the decrease of the 



27[. In most of non-relativistic 



Landau parameter Fr, to around —1 was rep orted in Ref. 

calculations ll^ ll^ ll^ l2l[ Q , decrease of the gap was reported. This 

is due to non-occurrence of liquid-gas instability and to the repulsive effect of the spin 
density fiuctuation. In the present model, the effect of the transverse uj polarization that 
corresponds to the spin density mode is weak. It is worth noting here that a recent ab initio 



calculation j23| concluded that the polarization effect in pure neutron matter is weak. 

It should be stressed that it is an approximation to consider only the medium (particle- 
hole) polarization in the relativistic model; the vacuum (particle-antiparticle) polarization 
should also be included. Since the gap is sensitive to the tiny remnant of the strong can- 
cellation mentioned above, future inclusion of the vacuum polarization would be important. 



This will be studied separately. According to Ref s. |42l l43|. the vacuum polarization leads to 
vector meson mass decrease. In a previous paper 4^ we examined this meson mass decrease 
using the in-medium Bonn potential and concluded that it reduces the gap. 

Other ingredients that are not considered in the present studvare inclusion of exchange 
of other mesons like vr and p, and the selfenergy effects . Study of superfluidity 

in the liquid-gas coexistent phase is also an interesting issue. 
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APPENDIX A: MEDIUM POLARIZATION INSERTIONS 



1. Feynman-density decomposition 
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2. Particle- hole-antiparticle decomposition 
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